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Abstract — Brain tumor segmentation for MR images is a 
difficult and challenging task due to variation in type, size, 
location and shape of tumors. This paper presents an 
efficient and fully automatic brain tumor segmentation 
technique. This proposed technique includes non local 
preprocessing, fuzzy intensification to enhance the quality 

of the MR images, k - means clustering method for brain 
tumor segmentation. The results are evaluated based on 
accuracy, sensitivity, specificity, false positive rate, false 
negative rate, Jaccard similarity metric and Dice 
coefficient. The preliminary results show 100% detection 
rate in all 20 test sets. 

Keywords- Brain Tumor, k-means clustering, Magnetic 
Resonance Imaging, Neutrosophic Set, Wiener. 

i. Introduction 

Magnetic resonance imaging (MRI) is widely used to visualize 
brain structures such as white matter, grey matter and 
ventricles, and to detect abnormalities in the brain. In general, 
different MRI sequences are acquired for better visualization 
of abnormalities with different contrast properties. Brain 
tumor segmentation in MRI is an important step in medical 
diagnosis because it provides information associated to 
anatomical structures as well as potential abnormal tissues 
necessary to treatment planning and patient follow-up [1]. 
Typically, radiologists perform manual tumor segmentation 
and impart the details to neurologists. This is a time- 
consuming process due to huge volume of images. In 
addition, the performance of manual segmentation is affected 
by the radiologists’ fatigue. Therefore, computer aided brain 
tumor segmentation from MRI data has been an active and 
promising research area in recent years. However, it is still a 
challenging problem to obtain the accurate and reproducible 
segmentation and characterization of tumors automatically due 
to the high variation in appearance of tumor tissues among 
different patients. 



Most of the existing techniques for brain tumor 
segmentation are not fully automatic due to expert-required 
model training or user-guided initialization or other user 
interaction, which leave significant space for increased 
automation, applicability and accuracy. The development of 
computer aided diagnosis (CAD) systems for brain tumor 
segmentation involves automatically assessing the probability 
of abnormality, and then highlighting the potentially abnormal 
regions for radiologists for further inspection. As a result, the 
probability of overlooking the tumor and the effort of manual 
examination of the whole 3D image are reduced. For this 
purpose, the CAD system should have the capability of fully 
automated segmentation of tumors. 

The existing work on automatic brain tumor segmentation 
methods are broadly classified into region-based and contour 
based methods. In the region-based methods, Clark et al 
proposed an automatic tumor segmentation from MRI using 
knowledge- based techniques with multispectral analysis [2]. 
First, initial segmentation is performed by an unsupervised 
clustering algorithm. The segmented image, along with cluster 
centers for each class is provided to a rule based expert system 
which extracts the intracranial region. Multispectral histogram 
analysis separates suspected tumor from the rest of the 
intracranial region, with region analysis used in performing 
final tumor labeling. Kaus et al proposed a method for 
automatic segmentation of small brain tumors using a 
statistical classification method and atlas registration [3]. 

Moon et al have used the expectation maximization (EM) 
algorithm and atlas prior information [4]. Prastawa et al 
proposed a brain tumor segmentation frame work based on 
EM algorithm and outlier detection, where tumors are 
considered as outliers of the Gaussian model. These 
techniques detect tumor regions using a registered brain atlas 
as a model for healthy brains [5]. Except for the professional 
work of atlas generation, large deformation of tumors requires 
significant modification of the atlas which may lead to poor 
results. Much more flexible techniques without priori 
information were developed to overcome these problems. 
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Philips et al first applied the fuzzy C-means clustering to 
segment brain tumors [6]. One of its drawbacks is that human 
interaction is required. Liu et al adapted the fuzzy 
connectedness framework for tumor segmentation with limited 
user interaction [7]. The fuzzy connectedness method has been 
applied by Moonis et al [8]. In this semi-automatic method, 
the user must select the region of the tumor. The calculation of 
the connectedness is achieved in this region and the tumor is 
delineated in 3D as a fuzzy connected object containing the 
seed points selected by the user. Dou et al proposed a fuzzy 
information fusion framework for brain tumor segmentation 
[9]. This method is sensitive to noise and needs some user 
interactions. Fully automated brain tumor segmentation from 
MR images based on improved fuzzy connectedness algorithm 
was proposed by Harati et al [10]. However, the fuzzy model 
works well only for hyper intensity (fully enhanced) tumors 
and exhibit poor performance on detecting non-enhanced 
tumors. 

In contour based methods, Lefohn et al applied level set 
implementation for brain tumor detection but for initialization 
of a first segmentation, in which user interaction is required 
[11]. Based on the visual inspection of the results of first 
segmentation, the level set parameters are tuned and the 
segmentation process is repeated. Ho et al incorporated region 
competition into level-set algorithm for brain tumor 
segmentation to overcome the initialization and weak edge 
leakage problems [12]. This method has high computational 
cost. Contour -based deformable models suffer from the 
difficulty of determining the initial contour, tuning the 
parameters and leakage in ill-defined edges. 

In this paper, a new method of automatic brain tumor 
segmentation technique is proposed based on neutrosophic K- 
means clustering method. First, the non local neutrosophic 
wiener filtering is applied to increase the quality of the image. 
Then fuzzy image enhancement technique is used to enhance 

the filtered image. Then finally k - means clustering method 
is used for brain tumor segmentation. 



ii. Proposed Method 

A. Preprocessing Using Neutrosophic Set 

Neutrosophy, a branch of philosophy introduced by 
Smarandache, as a generalization of dialectics, studies the 
origin, nature and scope of neutralities, as well as their 
interactions with different ideational spectra. In neutrosophy 
theory, every event has not only a certain degree of the truth, 
but also a falsity degree and an indeterminacy degree that have 
to be considered independently from each other [13]. Thus, a 

theory, event, concept, or entity, {S } is considered with its 
opposite [Anti-S] and the neutrality {Neut - S} . {Neut-S} is 
neither {s}nor{H«?;'-S}. The [Neut-S] and {Anti-S} are 
referred to as {./Vom -S} . According to this theory, every 
idea {H j is neutralized and balanced by {Anti -S} and 



{Non -5} [13]. NS provides a powerful tool to deal with 
indeterminacy. 



B. Non local Mean Algorithm 

The non local means (NLM) algorithm is applied to the 
noisy MRI in order to create the reference image. In the 
nonlocal means [14], given a discrete noisy image 
u - \u{i)\ i e /}, the estimated value 7VL [«](/) , for a pixel i , 

is computed as a weighted average of all the pixels in the 
image, 

NL[u](i) = Y MU j ) u(j) ( 1 ) 

jel 

where the family of weights {w(i,j)}j depend on the 

similarity between the pixels i and j, and satisfy the usual 
conditions 0 < w(i, j) < 1 and Yj Mi , ./) = 1 • 

The similarity between two pixels i and j depends on 
the similarity of the intensity gray level vectors w(AL)and 



u ( Nj ) , where N k denotes a square neighborhood of fixed 



size and centered at a pixel k. This similarity is measured as a 
decreasing function of the weighted Euclidean distance, 
II II 2 n 

w (N:) — u(N ■) , where Ct > U is the standard deviation of 

II J II 2 ,a 

the Gaussian kernel. The application of the Euclidean distance 
to the noisy neighborhoods raises the following equality 

EUN^-uiNjf =\\u(N i )-u(N j )f +2 a 2 

II J 112, a II J 112, a 

( 2 ) 

This equality shows the robustness of the NLM algorithm 
since in expectation the Euclidean distance conserves the 
order of similarity between pixels. 

The pixels with a similar grey level neighborhood to 

ll(N i ) have larger weights in the average. These weights are 



defined as, 



Mi, j) = 






z(0 



where Z(7) is the normalizing constant 



(3) 



\\u(N,-)-u(N j )\\ 2 a 

Z(i) = Y e h2 

j 



(4) 



The parameter h acts as a degree of filtering. It controls the 
decay of the exponential function and therefore the decay of 
the weights as a function of the Euclidean distances. 



C. The image in neutrosophic set 

The NLM filtered image is transformed into NS domain. A 
neutrosophic set image P NS is characterized by three 
membership sets T,I ,F . A pixel P in the image is described 
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as P(T, I, F ) and belongs to W in the following way: It is 
t true in the set, i indeterminate in the set, and f false in the 
set, where t varies in T , i varies in I and f varies in F . Then 
the pixel P(i,j ) in the image domain is transformed into the 
neutrosophic set 

domain P NS (i, j) = {T (i, j), I(i, j), F(i, j)}.T (i, j ), I(i, j) 
and F(i,j) are the probabilities belong to white pixels set, 
indeterminate set and non white pixels set respectively [15], 
which are defined as: 



T(iJ) 



g(zj)-gmin 
S max — S min 



(5) 



s{hj) = 
I(iJ) = 



i i+w/2 j+w/2 

-V" I X 

^ ^ m=i-w/2 n= j-w/2 

■/) ~ ^min 

c* c* 

^max ^min 



g(m,n) 



(6) 

(7) 



S(i, j ) = abs(g(i, j ) - g(i, j)) (8) 

F (hj) = \-T{i,j) (9) 

where g(i,j) is the local mean value of the pixels of the 
window. S(i, j) is the absolute value of difference between 
intensity g(i, j) and its local mean value g(i, j). 

D. Neutrosophic image entropy 

For an image, the entropy is utilized to evaluate the 
distribution of the gray levels. If the entropy is the maximum, 
the intensities might have equal probability. If the entropy is 
small, the intensity distribution is non-uniform. 

Neutrosophic entropy of an image is defined as the summation 
of the entropies of three subsets T ,1 and F [15]: 

En NS = En r +Erij +En F (10) 

max{r} 

En T =- 'YjPria) In p T (ci) (li) 

a=min{r} 

max{/} 

En, =- 2>/0) ln ^/(«) (12) 

<2=min{/} 

max{F} 

En p = - YjPf^) ^P F ( a ) (!3) 

a=min{F} 

where En T , Etij and En F are the entropies of sets 
T,I and ^respectively. p T (d) , p t {a) and p F (a) are 
the probabilities of elements in T, 1 and F respectively. 



E. a> -wiener filtering operation 

The values of / (i, j) is employed to measure the 
indeterminate degree of element P NS (i, j) . The set 
/ correlated with T and F , and the changes in T and F 
influence the distribution of element in 1 and the entropy of / . 

A co- wiener filtering operation for P NS , P NS [co) , is defined 
as [15]: 



P NS {co) = P(T(co), l{co), F(co)) 


(14) 


—• / \ f T I <co 


(15) 


Ta> J ) = wiener {T(m,n)} 
(m,n)eSjj 


(16) 


- , x f F I <co 

F(mh k 


(17) 


Fa 0'> J ) = wiener ( F(m , n)} 

{m,n)eS uj 


(18) 


j (j ; ) _ ^ T ^ f min 

1 CO j ) r» c 

Of ~ Of . 

T max T min 


(19) 


Sf O', j) = abs(T ( i , j ) - T ( i , j )) 


(20) 


— i i+w/2 j+w/2 

T{i,j)= X X T{m,n) 


(21) 



WXW m= i- w /2 n = j-w/2 



where 8j{i,j) is the absolute value of difference between 

intensity T(i,j) and its local mean value T(i,j) at 
(i, j ) after co - wiener filtering operation. 

F. Image enhancement 

Image enhancement is a basic method to improve the 
interpretability or perception of information in images for 
human viewers and it can provide better input for other 
automated image processing system. Here Fuzzy image 
enhancement technique is used for enhancing the tumor in 
MRI which is based on gray level mapping into fuzzy plane, 
using a membership transform function. The aim is to generate 
an image of higher contrast than the original one by giving the 
larger weight to the gray levels that are closer to the mean 
gray level of the image than to those that are farther from the 
mean. The intensification operation for the membership 
function is defined as [16], 

|2 p 2 (i,j) ju(i,j)<0.5 (22) 

[l — 2(1 -ju(i,j)) 2 ju(i,j)> 0.5 
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G. Segmentation 

For brain tumor segmentation, k -means clustering 
technique is used [17]. Clustering is the process of partitioning 
or grouping a given set of patterns into disjoint clusters. This 
is done such that patterns in the same cluster are similar and 
patterns belonging to two different clusters are different. For 

k -means clustering, the number of clusters k assumed to be 
fixed. Let the k prototypes (Wj , W, ,....W k ) be initialized to 
one of the n input patterns (Xj , X 2 ,....X n ) . Therefore, 

w i ~ x m -> / e {lj2,. ..k}, m e {l,2, ...n}. C , is the 

/^'cluster whose value is a disjoint subset of input patterns. 
The quality of the clustering is determined by the following 
error function. 

E ='t 

/=lx m eC, 



TN 

Specificity = 

TN + FP 



FP 

False positive rate (FNR) = 7 7 7 ; 

|77V| + |FP| 

= 1 -Specificity 

\FN\ 

False negative rate (FPR) = 7 7 — 7 ; 

\tp\ + \FN\ 

=1-Sensitivity 
Jaccard similarity metric 

\TP\ 

J(T,S) = 1 — 1 

TP + FP + FN 



Dice coefficient 



(25) 



(26) 



(27) 



(28) 



hi. Results and Discussion 

A. Materials 

The experimental datasets consist of 20 brain tumor 
patients MR images of T1 weighted and T2 weighted with 
different MRI pulse sequences. These images are obtained 
from PSG Institute of Medical Sciences and Research (PSG 
IMS & R), Coimbatore, Tamilnadu, India. Expert’s manual 
delineation of tumor in each MR image is treated as the 
ground truth. The method was implemented using MATLAB 
2010a (The Math Works, Inc) on the workstation with 
Pentium Dual core 2.4 GHz with 4GB of RAM and Windows 
XP 64-bit operating system. The segmentation results of the 
proposed method are given in Figure 1. 

B. Validation Strategies 

Several metrics are used for evaluation of the 
segmentations of same structure in an image [18]. Let X be 
the set of all voxels in the image. The ground truth is defined 
as T £ X , the set of voxels that were labeled as tumor by the 
expert. Similarly, S £ X as the set of voxels that were 
labeled as tumor by the proposed algorithm. The true positive 
set is defined asTP = T Pi S , i.e., the set of voxels 
common to T and S . The true negative is defined 

as TN = T O S , i.e., the set of voxels that were labeled as 
non-tumor in both sets. The false positive set is defined 

as FP = T O S , the non-tumor voxels are labeled as tumor 

voxels. The false negative set is defined as FN = T Pi S , 
the tumor voxels are labeled as non-tumor voxels. From these 
sets, success and error rates can be computed as 

\TP\ 

Sensitivity = 



TP 

D(T , S ) = 

— (JttpI + \fn\ + \tp\ + \fp\) 

2 (29) 

C. Results 

The detailed experimental results of Automatic Brain tumor 
segmentation of the proposed method is shown in Fig. 1. The 
Fig. 1(a) is the Original image and the filtered image using 
non local NS wiener filter is shown in Fig. 1(b). The enhanced 
image using Fuzzy intensification. Tumor traced by experts, 
and the Tumor segmented by the proposed method are shown 
in Fig. 1(c), (d) and (e) respectively. 

The evaluation result of the automatic segmentation of the 
proposed method is given in Table 1. 

iv. Conclusions 

In this paper, a fully automatic brain tumor segmentation 
technique is proposed based on neutrosophic preprocessing 

method and k -means clustering method. First, the non local 
neutrosophic wiener filtering is applied to increase the quality 
of the image. Then fuzzy image enhancement technique is 

used to enhance the filtered image. Then finally k - means 
clustering method is used for brain tumor segmentation. The 
preliminary results of the proposed show that the 100% 
detection rate in all 20 cases with average of high accuracy 
98.37%, high specificity 99.52% and lower missing rate 0.52. 



TP + FN 



(1A\ 
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(a) (b) (c) (d) (e) 



Figure 1 . Results of Automatic Brain tumor segmentation of the proposed method, (a) Original image, (b) Filtered image using non local NS wiener filter, (c) 
Enhanced image using Fuzzy intensification (d) Tumor traced by experts, (e) Tumor segmented by the proposed method. 
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TABLE I 

Results of Automatic Brain Tumor Segmentation of the Proposed Method 



Image 


Accuracy 


Sensitivity 


Specificity 


FPR 


FNR 


Jaccard 


Dice 
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0.8484 
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0.0118 
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0.7960 


0.8864 
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0.9874 
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0.9935 


0.0065 


0.1256 


0.8023 


0.8703 
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0.8957 


Case 12 


0.9784 


0.9659 
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0.0130 


0.0846 
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0.8864 
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0.9935 


0.0065 


0.1256 


0.8023 


0.8703 
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0.9956 
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0.9943 


0.0057 


0.0236 


0.9435 


0.9710 


Average 
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0.9926 


0.0052 


0.0985 


0.8446 


0.8981 
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